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Abstract 

We analyze the tubular phase of self-avoiding anisotropic crys- 
talline membranes. A careful analysis using renormalization group 
arguments together with symmetry requirements motivates the sim- 
plest form of the large-distance free energy describing fluctuations of 
tubular configurations. The non-self- avoiding limit of the model is 
shown to be exactly solvable. For the full self-avoiding model we 
compute the critical exponents using an e-expansion about the upper 
critical embedding dimension for general internal dimension D and 
embedding dimension d. We then exhibit various methods for reliably 
extrapolating to the physical point (D = 2, d = 3). Our most accurate 
estimates are v = 0.62 for the Flory exponent and £ = 0.80 for the 
roughness exponent. 
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1 Introduction 



The statistical mechanics of isotropic crystalline membranes has been the 
subject of much work in the last ten years [p], |2f. In the absence of self- 
avoidance there is a finite temperature crumpling transition from a low- 
temperature flat (orientationally-ordered) phase to a high-temperature crum- 
pled phase. The novel flat phase of phantom crystalline membranes is by now 
quite well understood, both qualitatively and quantitatively. The effect of 
self-avoidance on the phase diagram presents a much greater analytical and 
numerical challenge. While there is still some controversy, the bulk of evi- 
dence at present indicates that the crumpled phase disappears. It is possible, 
however, that this is the result of bending rigidity induced by next-to-nearest- 
neighbor excluded volume interactions. 

Rather surprisingly, it has been shown || that anisotropy has a remark- 
able effect on the global phase diagram of this class of membranes. For 
phantom membranes the flat and crumpled phases are isomorphic to those 
of the isotropic system (anisotropy is irrelevant in these phases) but there 
are intermediate tubular phases in which the membrane is ordered in one ex- 
tended direction (y) and crumpled in the remaining transverse directions (_L). 
Since self-avoidance is less constraining for configurations that are crumpled 
in one direction only, it is very likely that the tubular phase will survive in 
the more physical self- avoiding case, in contrast to the situation for isotropic 
membranes. Besides their intrinsic novelty, the study of membranes of this 
class may have important experimental and practical applications. First 
of all polymerized membranes with in-plane tilt order would have intrinsic 
anisotropy. In addition, polymerization in the presence of an applied electric 
field should produce anisotropic membranes H. 

The key critical exponents characterizing the tubular phase are the size 
(or Flory) exponent u, giving the scaling of the tubular diameter R g with the 
extended L y and transverse L± sizes of the membrane, and the roughness 
exponent ( associated with the growth of height fluctuations h rms (see Fig. [I]): 



Here Sr and Sh are scaling functions [Q, |5[] and z = K is the anisotropy 
exponent. In the phantom tubular phase (PTP) v and ( were computed in 
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Figure 1: A schematic illustration of a tubular configuration indicating the 
radius of gyration R g and the height fluctuations h rms . 

H, together with a self-consistent determination of the anomalous elasticity. 
The existence of the tubular phase has also been confirmed by numerical 
simulations |7| and the critical exponents measured are in excellent agreement 
with the theoretical predictions. In this paper we show that a careful analysis 
of the relevant operators in the free energy allows an a priori exact calculation 
of the anomalous elasticity as well as the above critical exponents. 

For self-avoiding membranes the model is much more difficult to treat 
analytically. By adapting the Edwards model for self-avoiding membranes 
to the geometry of the tubular phase, Radzihovsky and Toner H obtained a 
model free energy to describe this system. This was further studied by Bow- 
ick and Guitter ||, who utilized the multi- local-operator-product-expansion 
(MOPE) 0, 0] to perform an e = (d s c A - d) expansion about the upper crit- 
ical (embedding) dimension d^ A = 11. The phase diagram implied by this 
analysis is shown in Fig. Note the infrared stable fixed point (SAFP) with 
non-vanishing self-avoidance coupling b associated with the tubular phase. 
Bowick and Guitter also showed that the bending rigidity is not renormal- 
ized and computed the critical exponents to first order in e. They noted, 
however, that the extrapolation of these predictions to the physical tubule 
was not very robust against higher-order perturbations. 

Radzihovsky and Toner || have argued that the phase diagram described 
above is actually more complicated (see Fig. ||) for embedding dimension d 
less than a critical value d*, with > 3. They argue that the physics below 
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Figure 2: The phase diagram for self-avoiding anisotropic membranes with 
the Gaussian fixed point (GFP), the tubular phase fixed point (TPFP) and 
the self-avoidance fixed point (SAFP). 

d* is controlled by a new fixed point (BRFP) which is non-perturbative in e. 
This postulated fixed point is quite distinct physically from the SAFP. In par- 
ticular the bending rigidity picks up a non-zero anomalous exponent. Calcu- 
lating critical exponents at the putative BRFP would present the formidable 
challenge of a complete treatment of both self-avoidance and full non-linear 
elasticity. Reasonable estimates of v may, however, be obtained within the 
Flory approximation. 

In the present paper we begin with a careful analysis of the rotational 
symmetries of the tubular problem and their realization within a Wilsonian 
renormalization group approach M. This constrains the possible operators 
that may appear in the free energy and allows us to identify some operators 
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Figure 3: The phase diagram for self-avoiding anisotropic membranes with 
the Gaussian fixed point (GFP), the tubular phase fixed point (TPFP), 
the self-avoidance fixed point (SAFP) and the bending rigidity fixed point 
(BRFP). 

as definitely being irrelevant with respect to a broad category of fixed points. 
As a result of our analysis we can motivate the phase diagram Fig. |2|, which 
follows from the free energy studied in |J , with the incorporation of a relevant 
operator involving in-plane phonon excitations. The analysis of || assumes 
that non-linear elasticity terms are always irrelevant. It may therefore break 
down if new terms in the free energy alter the renormalization group flows. 
While this may change the character of the fixed point above, our analysis 
suggests that it is imperative to understand the SAFP in as much detail as 
possible. This is the focus of the present paper. Given the model we next turn 
to the actual calculation of reliable critical exponents in the tubular phase. 
This is done by generalizing the calculation of || to manifolds of arbitrary 
internal dimension D embedded in general dimension d. We analyze a class 
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of generalized e-expansions that allow us to determine an optimal path from 
the line e = to the physical point (D = 2, d — 3). Our most accurate 
estimates are 

v = 0.62 

C = 0.80. (2) 

Furthermore, we show that the critical exponents determined in this method 
are extremely close to the Flory prediction, particularly for d > 3. This may 
be regarded as strengthening the predictions of the otherwise uncontrolled 
Flory approximation. 

The outline of our paper is as follows. The model is described in Sec. |2] 
along with an analysis of its symmetries and their implementation in a Wilso- 
nian renormalization group framework. This leads to a clarification of the 
global phase diagram and a proposal for the simplest free energy capturing 
the essential large distance physics of the tubular phase. This is followed 
in Sec. [5] by a derivation of the scaling relations connecting the fundamental 
critical exponents. The special case of the phantom tubule is treated in detail 
in Sec. |j. The full physical problem of the self-avoiding tubule is tackled in 
Sec. H|, where critical exponents are computed via a generalized e-expansion. 
We also compute corrections to the Flory and Gaussian variational approx- 
imations. A brief summary of our results is given in Sec. |6|. Finally, some 
technical details of the e-expansion are left to the Appendix Sec. [A]. 

2 Model 

A membrane configuration may be characterized by giving the position r(x), 
in the <i-dimensional embedding space, of a point in the membrane labeled 
by a .D-dimensional internal coordinate x. A physical membrane corresponds 
to the case d = 3 and D = 2. 

In |3|, Q the most general Landau- Ginzburg- Wilson free energy F for this 
system is constructed by expanding F to leading order in powers of r(x) 
and its gradients with respect to internal space x, taking into account global 
translation and rotational invariance. We will consider the case in which the 
membrane is isotropic in D — 1 membrane directions (denoted xj_) orthogonal 
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to a distinguished direction y. The resultant free energy is given by 

1 



F(f(x)) 



x±dy n ± (d ± r) + K y (d r) 



+K ±y d 2 y f- d 2 ± r + t ± (d^f) 2 + t y (d y f) 2 



+ ^r(d^r. + -f(d y f- d y r) 

+u ±y (d^f- d y f) 2 + ^(^f • dir) 

+ ^Jd D xJ d D x'5 d (f(x) - f(x')), 



(3) 



where the parameters denote bending and elastic moduli. Note the complex- 
ity of this model - it has eleven free parameters. In mean field theory the 
non-self-avoiding limit (b=0) yields a phase diagram with flat and crumpled 
phases separated by a tubular phase 0. 

In this paper we will be mainly concerned with the tubular phase (TP) 
beyond mean field theory. In this case we may expand f in the Monge 
representation: 

f(x) = (C^ + n(x),Mx)). (4) 

The free energy is now a function of u and h. Before simplifying Eq. [5] let us 
discuss the symmetries of the tubular phase. 

Since the free energy must be invariant under global rotations of the 
tubule it is expressible in terms of the complete set of tubular rotationally 
invariant operators. These are 



E(u,h) 

F a (u,h) 
F a p(u, h) 
G y (u, h) 
G ya /3(u, h) 

Indeed, Eq. [5] becomes 

1 

2 



F{u,h) 



1 1 

d v u + 2i d v h ) 2 + 2^ M ) 2 

d a u + d y hd a h + d y ud a u 
daudpu + dahdph 

(d 2 yU) 2 + (d 2 h) 2 

(d 2 U)(d a pU) + d 2 y hd a ph . 



(1 X±dl/ 2(y(ty + Uyy(y)E(U 1 fl) + Ky( y Gy(U, fl) 



(5) 
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+K yl _G« a (u, h) + 2u yy ^E 2 (u, h) + {t ± + v ±v )QF2 + C A y u Xy F a F a 

+2v ±y e y E( U , h)F: + ^FfFS + v -f{F:f 

+ b - J d D ~Wdyd D - l ^ ± dy'5 d -\h(-K L ,y) - h(x' ± ,y')) 

x$((y(y-y r ) +u(x±,y) -u(x.' ± ,y')) (6) 

Since we are interested in the critical properties of the free energy Eq. ^|we 
may simplify by dropping irrelevant terms. Simple power counting around 
the Gaussian fixed point is usually enough to determine the relevancy of 
operators but in this case the situation is more involved and requires a careful 
analysis of the symmetries of the problem, to which we turn now. 



2.1 Wilson RG in the tubular phase 

We apply the Renormalization Group (RG) a la Wilson || to the free energy 
Eq. ^| While this approach is usually more involved for extracting actual 
numbers than the more conventional field theory approach [IU|, it is more 
general and allows an easier analysis of the irrelevant operators, key to de- 
ciding which terms to retain in the free energy. The crucial point in Wilson 
RG is the RG transformation. This is a two step procedure: the blocking and 
the rescaling. 

There is considerable freedom in the choice of blocking. We chose deci- 
mation in momentum space, where in order to simplify the calculations an 
anisotropic spherical momentum regularization is assumed. The blocking just 
consists in integrating over an anisotropic shell of thickness e~ l , I G [0, +oo). 
That is, 

e -F l (u,h) = f jj du( q± ,q y )dh(q ± ,q y )e- F ^ . (7) 
{|gx Utiles 

The region B consists of three sectors 

!1 > |g_J > e~ l , e~ zl > \q y \ > 
1 > \q±\ > e~ l , 1 > \q y \ > e~ zl (8) 
e~ l < \q±\ < , 1 > \q y \ > e~ zl 

where the exponent z accounts for the anisotropy of the system. This block- 
ing is very similar to the one used in 0. 
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The rescaling is anisotropic as well and is given by, 

q' ± = e l q ± h'{q') = e^ D - 1+z+ ^ l h{q) 

q' y = e zl q y u'(q') = e^ D - 1+2 ^ l u(q) , [ ' 

where v is the other exponent that appears in the theory. 

The result of performing a renormalization group transformation up to 
time T, is the Wilsonian free energy 

Fi(u',h') , (10) 

where the u' and h' fields have the same range as the original ones. The free 
energy evaluated at / = is, by definition, Eq. §. 

For future reference, let us work out the simplest fixed point in Eq. ^|, 
the Gaussian fixed point. Although this fixed point is not of direct physical 
interest it plays a central role in many considerations (see Fig. ||). This 
fixed point may be studied by retaining only the quadratic terms in the free 
energy Eq. |^, and applying the RG transformation just defined. We easily 
get (hereafter dropping all primes in the rescaling) 
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+ (e^- 1 - 3 ^ 4 ^ + e^- z+ ^ l g ± ql)u(-q)u(q)} . (11) 

Imposing that the Gaussian fixed point is given by the terms involving h, 
the exponents z and v are readily computed to be 

z = - , 2v=--D (12) 

and the exponents for the operators associated with the couplings are uniquely 
determined. The Gaussian fixed point is thus g y — g± — 0. The coupling g± 
defines an irrelevant direction for D > 3/2, with exponent | — D, while g y 
defines a relevant direction for D < 5/2, with exponent \ — D. The Gaussian 
fixed point is therefore infrared unstable. 

2.2 The rotations of the tubule 

For the general free energy of Eq. |6| the rotations of the tubule are imple- 
mented by 

u — > u cos 6 -\- sin Oh -\- (cos 6 — 1)?/ 
h — > h cos 6 — sin 6u — sin Oy 
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where we have simplified by rotating just one component of h. The symmetry 
transformation above is unusual in that it changes under the action of the 
renormalization group. This happens because rotations of the tubule mix 
two sets of fields - the in-plane and out-of-plane phonons - having different 
scaling dimensions. In fact, it is straightforward to show that Eq. |H| is 
realized at time T by 

u — > ucos9 + e-^- z Hsme + e' 2 ^-^ l (cos9-l)y 

h cos 9 -e^- z H sin 9 -e'^ 1 sin 9y . ^ ' 



The above transformation is an exact symmetry of the free energy Eq. [10 
This transformation depends explicitly on I and prevents a simple construc- 
tion of invariant free energies. At large I, however, we may derive an l- 
independent version. Define 9 = Ae^~^ 1 and assume that the condition 

u(l) - z(l) < (15) 

is satisfied. Near the fixed point, scaling relations to be derived later show 
that 

v-z = -(v-D + l) (16) 
3 

and therefore v — z < for all v < D — 1. The physical case D = 2 requires 
v < 1, which is always valid. 
Eq. [14| is then, for large /, 

u — ► u + Ah - \A 2 y + 0(e 2 ^ 1 ) 

h -> h- Ay + 0{e 2 ^ 1 ) . [ ' 

The generalization of this symmetry to an arbitrary rotation involving h is 

u + Ah - \A 2 y + 0{e 2 ^- z ») , , 

h^ ti-Ay + 0(e 2 ^ 1 ) , 1 ' 

which is the tubular phase version of a symmetry noted earlier in || for the 
free energy describing the large distance properties of the flat phase. 

2.3 The large distance free energy of the Phantom 
tubule 

Let us apply the previous considerations to the construction of the free energy 
for the large distance properties of phantom tubules (Eq. ^| with 6 = 0). 
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In [[| ||, the free energy 



F(u, h) 



- / d°-Wdy K{d 2 v h) 2 + t{d a hf 



+g±{d a u) 2 + g y {d y u + -(d y h) 2 ) 2 



(19) 



is given as that describing the right large distance properties of the TP. 

The first thing to notice is that this free energy is not invariant under the 
symmetry Eq. [L8|. The free energy with the correct invariances is given by 

F(u,h) = - [ d D -Wdy\n(d 2 h) 2 + t{d a hf 



(20) 



+gj_{d a u + d a hd y h) 2 
+ 9y (d y u + ~(d y h) 2 ) 2 



since the operator d a u + d a hd y h is rotationally invariant. 

It is important at this point to recall that the symmetry Eq. [18] is exact 
up to 'irrelevant' terms, and the coupling gj_ is irrelevant for all the entire 
range of D (including D = 2) in which the TP exists. If we therefore insist 
on including irrelevant operators around the Gaussian fixed point, our free 
energy would certainly contain a non-invariant term under Eq. [TR 



F(u, h) 



- I d D ^ ± dy \K(d 2 h) 2 + t{d a hf 



+g i \_\d a u + dahdyh) 2 + g [ ^'(d a u 



(2), 



+ g y (d y u + -(dyh) 



2\2 



(21) 



Indeed, this is the combination that appears, up to higher irrelevant terms, 
in the general expression for the free energy Eq. || as g^ is the coupling to 

(2) 

the F a F a operator, and g ± is the coupling to F£. 

The usual strategy, nevertheless, is to keep just those operators that de- 
fine relevant directions of the Gaussian fixed point. It is these directions 
that flow towards new infrared fixed points, unless a first order transition 
occurs. Adopting this approach the relevant free energy for the phantom 
tubule would be 



F(u, h) 



d u -^ ± dy K(d 2 h) 2 + t(d Q hY 
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+ 9«(d y u + ~(d y h) 2 ) 2 



(22) 



where g y defines a relevant direction for D < 5/2 which terminates in the 
tubular phase fixed point (TPFP) as shown in Fig. |2|. Note that the symmetry 
Eq. O is, indeed, preserved. 



2.4 The large distance free energy for the self- avoiding 
tubule 



Now let us return to the more physical model with the self-avoidance term, 
J dyd D ^ ± J dy'd - 1 ^^, y) - f(x' ± , y')). (23) 



b 
2 

restored. 



Following the discussion in subsection 2.2, we simplify the self- avoiding 



term Eq. £3] by demanding invariance under the symmetry Eq. [18], 

b - J dydP-hcidP-V^ih^y) - h(*' x ,y)) , (24) 

with irrelevant terms dropped. 

The scaling dimension of the new perturbation Eq. ^3] at the Gaussian 
fixed point is e = 3D — | — (| — D)d. The b coupling therefore defines a new 
relevant direction for D tubules embedded in dimensions d < d^ A , where 

d SA {D) = QD-^ 
c K ' 5-2D v ' 

Below the upper critical dimension d^ A the Gaussian fixed point is infrared 
unstable under this perturbation, and the large distance properties of the self- 
avoiding tubule are described by a new fixed point (SAFP). This new fixed 
point merges with the Gaussian fixed point at the upper critical dimension 
where self-avoidance becomes a marginal perturbation. We therefore expect 
the critical properties of the self-avoiding tubule to be perturbative in e, as 
pointed out first in (see Fig. 

In 0, however, it is claimed that this simple scenario is valid only for 
tubules embedded in dimensions d close to d^ A (D). For any dimension d 
lower than d* (where <i* < A ), they argue for the existence of a distinct fixed 
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point, the bending rigidity fixed point (BRFP) (see Fig. |3[). This fixed point 
is postulated to describe the actual critical properties of the self-avoiding 
tubule for the regime d < d„, including the physical case of the D = 2 tubule 
embedded in d = 3. If this scenario is true, the critical properties of the 
self-avoiding tubule are not perturbative in e. Analytical predictions become 
then extremely difficult, as there is no evident small perturbative parameter. 

At this stage, therefore, we need to understand better the topology of 
the RG flows in the case where self-avoidance is included. Let us review 
the arguments of [f|. They consider the free energy Eq. [19|, together with 
the self-avoiding term Eq. They include all relevant directions from the 
Gaussian fixed point, and an irrelevant one defined by g±. They apply the 
infinitesimal renormalization group a la Wilson to derive an equation for the 
evolution of couplings. The crucial equation in their analysis is the RG flow 
equation for g± 

^=[Au-z + D-3]g ± . (26) 

Now, as the RG is iterated starting near the Gaussian fixed point, g± de- 
creases to zero while the rescalings u(l) and z(l) flow towards their SAFP 
values. For sufficiently small embedding dimension d and large enough I the 
sign of the /3-function for g± changes sign. The coupling g± then flows to the 
BRFP g* ± = oo (see Fig. |3|). This argument can be made more quantitative. 
Under very reasonable assumptions, Eq. ^6] leads to a lower bound for d*, 
the highest embedding dimension in which the BRFP prevails, 

AD - 1 

d*{D) > T3 -p • (27) 

In particular, d*(2) > 7/2 > 3, so the physical tubule (D = 2, d = 3) is, 
according to ||], described by the BRFP. 
It is apparent that the operator 

d a ud a u (28) 

plays a fundamental role in this argument. Let us examine it more closely. 
In an expansion in irrelevant operators around the Gaussian fixed point, it 
appears in two ways, which we labeled g± and g± in Eq. |2l]. 

First of all, the operator associated to g± is invariant under the symmetry 
Eq. as it appears in the invariant combination 

d a u + d a hdyh . (29) 
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(2) 

In contrast g ± couples to a subdominant piece of the operator F£ = d a ud a u+ 
d a hd a h (see Eq. [5]). In fact, from our earlier symmetry arguments, it is sup- 
pressed by a factor 0(e 2 ^ u ~ z ^) with respect to the dominant piece (d a h) 2 
which couples to the marginal direction t. Provided v — z < the coupling 
g ± is thus irrelevant and may be dropped from the free energy. 

We have argued that the most general free energy dictating the large dis- 
tance properties of the tubule is given by Eq. ^ together with self-avoidance 
(Eq. |24j). For g± vanishing, the infrared stable fixed point of the theory is the 
SAFP. The key issue is now whether this fixed point is stable with respect to 
perturbations by g±. Since the properties of the SAFP are perturbative in 
e, the same applies to the critical exponents. Experience with typical multi- 
critical behavior suggest that we should not expect the exponent associated 
with the g± direction to change so much from its gaussian value 3/2 — D that 
it changes sign [[K] . 



In conclusion, the simplest free energy describing the large distance prop- 
erties of the self-avoiding tubule is given by 



K{d 2 y hf + t(d a h) 2 + 9y (d y U + lidyhf" 1 



F(u,h) = ^Jd D -^ ± dy 

+ h - J dyd^Wd - 1 ^- 1 ^^) - h(x! ± ,y)) ■ (30) 
This is the starting point of all our subsequent analysis. 



3 The Scaling Relations 

Having identified the right free energy, we turn now to the derivation of 
the different critical exponents of the theory. We use the conventional field 
theory formalism [1C], following 0. 

The scaling dimensions of the fields and coordinates are [y] = l,[x±] = 2, 
[h] = \-D and [u] = 4 - 2D. This implies 

[b] = -e , [g y ] =2D-5, (31) 

with 

1 5 - 2D , 
£ = 3D ~ 2 2~ d - (32) 
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Following the arguments in Jj], one can show that the free energy Eq. |30| 
renormalizes onto itself with 



F(u, h) = ~ / dP-^ldy \ZK(d 2 v h R ) 2 + Z ± t(d R h R ) 2 



b R Z b fi £ 



(dyU« + ~(d y h R ff 



(33) 



d D - 1 ^d D - 1 ^dy6 d - 1 (h R (^,y) - h H (xf,y)) 



R(-JR 



where the Ward identity implied by Eq. [IB] is used so that there is no inde- 
pendent wave function renormalization for the field u. Furthermore, it is not 



difficult to show, using the MOPE formalism Hl2|| , that the bending rigidity 
is not renormalized so Z — 1, as first pointed out in 0. Thus we have 



h R (x R , y) 

b R 



Z ± 4 h{x ± ,y) 

zl* ± 
b^z^z^- D ^ 



(34) 



g- = fi 5 ~ 2D g y Z 



D-l 
2 



Using these definitions we will consider two correlators, which enable us 
to determine the exponents of the theory. In the original paper ||, the 
correlator 

^(xjl, y) = - 27rf3T) ( (/I(x± ' y) - K(0 ' 0))2 ) (35) 

was considered as well the correlator involving the u fields, 
G u (x±,y) = (d y u(x.±,y)dyu(O,0)) . 
At the fixed point, the first correlator satisfies 

D-l 



g_ + g x 9 
dji 2 L <9x_i_ 



which, combined with simple scaling law 

d 



5}G H (x ± ,y) = 



d d 

^tt - yir - 2x± « — 

ojj, ay axj_ 



+ (5-2D)\G R (^,y) = 



(36) 



(37) 



(38) 
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gives us the fixed point renormalization group equation, 



4y + ¥^- 2( } G ^« ) = °- 



(39) 



A renormalization group equation may also be derived for G u . To do 
so we must use once again the Ward identity that fixes the wave function 
renormalization for u. Eq. K7] is now, 



Eq |38| reads for the u case, 



d d n d 
A*7j 2xj_— — 
on ay axx 



(6-4D)\G%{x ± ,y) = 



leading finally to 



where 



A 

2 



logZ± 



o 



c 



4 + 5 

5 - 2D 1 - D 



1 + 



2 4 
1 - D 



5 



z 



Both, Eqs |39| and |2 may be solved explicitly, yielding 

G h (x±,y) = 

G u (x±,y) = y 2Cu F[ { - t r ) = | x ± 



y K Fi(—^—_ 
l x ±r 
y 



X -L 



X -L 



x ± 



(40) 



(41) 



(42) 



(43) 



(44) 
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where v = * . Transforming Eq. ^ to momentum space gives 

G^q)- 1 = \px\ 2+ ^f(^-r) (45) 

q 



G u (p±,q)- 1 = \p±\ z ^h(-±-) 

P-L 



with 



?7_l = -2 + 4z 

% = -■ (46) 

z 

These scaling laws were first derived in || and || respectively. 

We conclude that all the critical exponents of our Free energy Eq. ^ at 
any putative fixed point may be expressed in terms of a single parameter, 
say 5. The task of computing critical exponents translates into the task of 
evaluating 5 at the corresponding fixed point. 



4 The Phantom tubule 



The theoretical considerations in subsection |2.3| lead us to consider Eq. |22| as 
the right free energy describing the large distance properties of the phantom 
tubule. In fact they allow us to solve the phantom tubule phase exactly, 
simply by performing the shift 

u^u'- \ ( V dz(d z h? , (47) 
I Jo 

where the lower bound for the integral is arbitrary and corresponds to trans- 
lations of the zero mode. The free energy Eq. ^ is then a sum of Gaussian 
terms. Let us compute the anomalous elasticity, determined by the correlator 
Eq. 

G M (x - z) = (d v u{x)d y u{z)) , G u (p) = J d(x - z)e- p ( x - z )G(x - z). (48) 
The elasticity constant is given by 
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At tree level g y (p) = g y - 

The general case amounts to performing the shift Eq. 47. Equivalently 
the loop expansion may be performed to all orders. The diagrams consists 
of a necklace of h loops which can be resummed, yielding 

1 1 , d-1 r q 2 y (p-q)l 

- — + — — / d q±dq y - 



g y (p) g y 2 J ^ + t<&)(K(q y -p v )4 + t(<i L -P±) 2 ) 



which for p 7^ is 



^(P) = Ta — ~\\ , (51) 

1 + &$tLf{D - \) V f - 5 C(^) 



where f(d)= J d d i 



(* 2 +i) 4 ' 

( D-5 
^^ + j(x(l-z) i +(l-x)z 4 )\ 2 , (52) 

and the exponent z is 

2 = 1/2. (53) 

Recall that Eq. |5l] is valid for any value of g y . The Gaussian fixed point 
(g y = 0) is unstable to perturbations along this direction, and the coupling 
g y is driven to g y — 00 in the infrared, which is the fixed point describing 
the physics of the phantom tubule (PTFP). At the PTFP g y (p) has the form 
P y w g(r^jr^), as predicted by Eq. (44] at any fixed point. 

We can easily recover now the results in || H from our exact solution 
Eq. and Eq. ^ at the PTFP: 

• For y = — > 00 



K,n- B 2 f+00 (1 _ z )2(D-3) _ z 2(D-3) 



Ipj-I 

we have 

(54) 

which converges, both in the infrared and the ultraviolet for | < D < | 
(the dimensions in which the tubular phase exists). For small p y 

9y(p)-P v i ,Vu = 5-2D. (55) 
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Py 



For y — | — 
we have 



-2D+5 



C(y) y ^ ~ V— — {-) D f (x(l-x)) 2J ^ [ °° dz^l + z)™ (56) 
2 k Jo Jo 

so C(y)^Q ~ y Vu Const. where Const, is a convergent integral for | < 
D < § . For small p ± at D = 2 

ft,(p)~pi /a - (57) 

To conclude let us connect with the results in section [| At the PTFP we 
have 6 = 0, and rest of the exponents from this result and scaling relations. 

5 The self-avoiding tubule 

We have argued that Eq. ^FJ is the appropriate free energy to consider once 
self-avoidance is included. The task of computing the critical exponents of 
this theory is not easy, since from we know that the results of the e(d)- 
expansion are not very robust to higher order perturbations. 

In order to get an estimate for the exponents at the SAFP, we compute the 
critical exponents to lowest order in e(D, d) for arbitrary internal dimension 
D. Existing techniques [13|, [14|] then allow us to perform more sophisticated 
extrapolations which produce reliable estimates for the critical exponents. 

5.1 The computation of 5 



We follow the MOPE formalism [Lj|, employing dimensional regularization 
and minimal subtraction, used first in this problem in [H], to compute the S 
exponent. 

Within the MOPE formalism, one may prove that the free energy Eq. |SD| 
renormalizes onto itself. It also identifies the diagrams to compute that yield 
the RG functions determining the critical exponents. For details on this 



formalism we refer to [11, 12], and for its implementation in the tubular 



case, we refer to the original BG || calculation. 
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The first step is to compute the two-point function Gh for arbitrary D at 
6 = 0. The result is 



\2-D 

G° h (x±,y) 



1 r+°° Di i 
|x_l|2 / dtt* Ki-d (t) cos(t 2 uO 
, 2 - j D)(2tt)^ L 1 X| Jo — W V ; 

+^ / dtt—K 3 -p(t) sm(t*w) (58) 

Z Jo 2 

where u> = — ^-r and is a modified Bessel function. There are two par- 

|xj.|2 

ticular cases of interest. At y = we have 

^(xx,0) = - |X± 'f D ^ ) ^: i) - (59) 
(I - D)tt— 25 

For the physical value -D = 2 it follows from Ki(t) = K_i(t) = (j^J e~*, 
that 

^(x,y) = -j^e- 2 #-ferfS, (60) 

27T2 4 I 

where erf(x) denotes the error function. This result is in complete agreement 
with that quoted in ||. The next step is to perform the MOPE for the 
operator 

<j>{^ ± ,y} = 5^ 1 \h^ ± ,y)-h(^ ± ,y)) . (61) 
This is easily done using standard techniques, with the result 

0{x ± ,x' ± ,y} = C*(x ± -x' ± )+Cf (x ± -x' ± ) : V^(x° ,y)V^(x° ,y) :+•-., 

(62) 

where x° ± = ± 2 x , and the Wilson coefficients are 

1 u a u^ 



(4vr)^(-G°(u,0))^ ' ^ 4(4 7 r) i ¥(-G'£(u,0))^ ' 

(63) 

We also need the MOPE for the product of two of these operators. One finds 

0{*r, z i, yi}H*2, 2/2} = c^(x^-j^, z^-z^, y 1 -y 2 )0{x ± , z ± , y}+- ■ ■ , 

(64) 
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where xj_ = Xl ^ X2 , y = yi + y2 and z± = Zl t, 7 " 2 and 

Ctju, V, w) 



1 1 



(47r)V(-G°(u, W )-G'0(v,«;))V 



(65) 



This is all we need to compute the critical exponents. This MOPE corre- 
sponds to the diagrams in Fig. [|. The last diagram for the renormalization 
of b is not necessary to compute, as it cancels against the renormalization of 
Z,. 




y 2 x 2<^ 




Z3-X 2 



Figure 4: One loop diagrams contributing to the renormalization of the free 
energy. 

Expanding the renormalized action Eq. |33], and using the MOPE Eq. 
b R n 



2 

b R /i £ 



J d D - 1 x ± d I3 - 1 x , ± d|/0{x ± ,xl,y} 
J d°- 1 x ± d D - 1 x / ±( iyCj(x ± - x' ± 



(66) 



+ J d B -Vd fl_1 xidyCj' /J (x J .-x , ± )V a K(x5.,y)V^(xi,y) + 



The first term of Eq. |66| renormalizes the identity operator and therefore may 
be neglected in computing expectation values of operators. The second term 
determines Z\, with the result, 



1 + 



1 



D — 1 2 



3D+3 11-3-D 
7J- 10— ID 2 5-2D 



6 R 1 



(r(t)r(f -l)) 5 - 2D 



20+2 -p/_D_i- 



(67) 
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At one loop there is no renormalization for Z, a result that is also true at 
any order in perturbation theory ||. 

Expanding the (^-function, and performing the MOPE in Eq. we find 

x J d^zd^wdyC^z^y) , (68) 

where higher terms in the MOPE Eq. ^ are neglected as they do not give rise 
to poles in e. To find Z b , we must compute the last integral in Eq. This 
is done by performing the angular integration and then changing variables 
to u = — — and v = — — . The result is 



d-l 



j d^zd^dyCUz^y) = -^rj^yrr) (69) 



x(^-D){2ir)^y j'dyy 
x / du dv 



J_ i 1( 2D-3 7 ,2D-3 



where 



(/(«) + /(«))■ 

t l/2 



-1 j 
2 



f( u ) = v?- 2D \u dtt^K^nU) cos( ) 

y Jo 2 u 

/ dtt^ATa^t) sin( H. (70) 

Jo 2 w J 



1 /-+oo 
+ - 



Using 



J_ „_!_ u 2D-3 v 2D-3 

iy du / ™ (iv t— j- (71) 

'o (/(«)+ /(t;))^ 



-x /-+oo u 2D-3 v 2D-3 



•A> (/(«) + 

where p(y) is a continuous function that vanishes at y — 0, and adding a 
factor of 2 corresponding to the two ways one can perform the MOPE in the 
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diagram in Fig. 0, we get finally 



, 4D-3 (4P-3)(P-3) . . 

(|_D)B=au2 2 < 5 - 2D ) + , 



-7(D-1) 

7r ^w r (^i)2 



with 



7(D) 



du 



dv- 



u 2D-3 v 2D-3 



(72) 



(73) 



(/(«) + /(«))™ 

We have thus succeeded in renormalizing the theory at the one loop level. 
The evaluation of this integral is discussed in the Appendix. The next step 
is to compute the exponent 5. We begin with the computation of the j3 



function. There are two of them. Defining a\ and b\ via Z_\_ = 1 + — a\ and 



1 + — &i, we have 



Mb 1 



-eb R + (ai + 



7-e 



and 



2(5 - 2D 
D - 1 



-h)(b 



R\2 



^(b R ,g R ) = (5 - 2D + r_^A(0)<? y . 



(74) 



(75) 



where A(6 ) = /j,j-logZ±. There is a nontrivial fixed point (the SAFP) at 



and for D < 5 - A (^*)/ 2 
ana, tor u < 2 _ A(bR ,y 2 



di + 



R 



7-e 



2(5-2L>) 



6i 



<7„ = +oo . 



Using 



5 = A(b R *) = fi^-logZ„ 



d 



dfx 



b R =b R * 



we obtain the final result 



< i_ »i 

2(5-2D) ' 6i 

Plugging Eq. |7| and Eq. |7| into Eq. |7| yields 



5 



f +tf(D)7(D) 



(76) 

(77) 
(78) 

(79) 
(80) 
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where 

2D+2 2Q 2 -25U/2+27/2 



m = p-i) - - 2 4 r %-i, ■ (si) 

nT{— ) 

The scaling relations after Eq. ^2] and Eq. ^6] determine the rest of the critical 
exponents provided we have a good determination of 5. As only the first term 
in the e expansion is available this will require refined methods to improve 
the perturbative expansion. This is an involved subject to which we now 
turn. 



5.2 Analysis of the results 

From Eq. ^ and the scaling relations we get explicit forms for the critical 
exponents of the self-avoiding tubular phase. For example, the radius of 
gyration or size exponent v reads, 

^ — 2D 

u(D) = + ui(D)e(D, d) + ■ ■ ■ , (82) 

where vi{D) is plotted function of D in Fig. [5| 

As already noticed in 0, a direct application of Eq. ^ to a physical 
membrane is not robust with respect to second order corrections. This is a 
consequence of the point (D — 2, d — 3) being too far from the point (2, 11) 
on the critical curve e = 0. One can try, instead, to perform a generalized e 
expansion around any other point on the critical curve (D , d = frf^) (see 
Fig. |6|) and hope to find a new expansion in which the corrections to Eq. |82| 
are minimized. In this case one may expect reliable one loop results. 



As an example, let us rewrite Eq. |82| in terms of D — D , and keep the 
leading terms, 

u(D) = b -^^- I ^^ + u 1 {D Q + D-D )E{D,d) + O{e 2 ) 

= 5 -^-^^ + MDo)e(D,d) (83) 
+ O(e 2 ,6(D-D ),(D-D ) 2 ) . 

One can expand around any point (D , do) in the e = curve, but at 
the expense of dealing with the double expansion in e and D in Eq. R31 
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Figure 5: Plot of Vi{D) as a function of D. 



Furthermore, critical quantities depend, in principle, on a new parame- 
ter D , but should obviously be independent of it. There are established 
techniques to select the best Do, such as the minimal sensitivity method 
of Hwa ||13|| . Anyway, the expansion in e, D is just a particular case of 



a more general situation | 14| , as we may choose any new set of variables 
{x(D,e),y(D,e)} and re-express Eq. ^3] as an expansion around the critical 
curve (x = x(D , 0),y = y(D , 0)), 



u(D, e) = v(x, y) = v(x , y ) + Axu lfi (x , y ) + Ayis 0A (x , y ) + 



(84) 



where Ax = x(D, e) — xq, Ay = y(D, e) — yo. 

Now the goal is to choose a good set of variables {x, y} from the infinite 
number of possibilities. 

A good expansion variable must meet two basic requirements; 
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Figure 6: The solid line is the e = curve. The naive e expansion is marked 
with a dotted line. Other possible expansions from the physical interesting 
case (D — 2, d — 3), are marked with dashed lines. 
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1. Independence of -Do, at least in some interval of Dqs. 



2. Within this region, results must be quantitatively correct. 

The way we choose to meet requirement 1, is to look for v to plateau 
within some given accuracy for an interval of Dq. To meet requirement 2, we 
will systematically study the (2, d) line, starting close to the d = 11 where 
even the naive e expansion is expected to work and agree well with the Flory 
estimate. At large d a good expansion variable should exhibit a broad plateau 
around the Flory result. At this stage, the expansions that do not deliver 
accurate enough results cannot be trusted, and should therefore be rejected. 
As d is decreased we check that the plateau remains stable for the good 
expansion variables, so that we have a reliable extrapolation at d — 3. In 
addition, we expect reasonable agreement with the Flory estimate. 

We examine the v exponent as well as the ( exponent using the above 
techniques. The reason for analyzing both exponents is that they have dif- 
ferent dependencies on 5, as expressed in Eq. [K| This may result in different 
expansion variables being appropriate for different exponents. Of course, 
results must be consistent with the scaling relations expressed by Eq. [|3| 

In this paper, we implement Hwa minimal sensitivity scheme ||13|| , and we 
explore the following distinct expansions; 



{x 


= D,y 


= e}. 


{x 


= D,y 


= d}. 


{x 


= D,y 


= D (d) 


{x 


= £,V = 


= D (d)} 



5d+l 
2(3+d) 



}• 



which have been previously used in a different context in 14 . 



5.2.1 Corrections to mean field 

The analysis leading to the extrapolation for the v exponent, may be sum- 
marized in Fig. 0. At (D = 2, d — 8) all sets of variables A,B,C,D give 
consistent results. Nevertheless, the expansion D shows the flattest plateau 
which is in complete agreement with the Flory estimate. This singles out 
expansion D as the best, and in fact we have taken as the actual v its value 
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at the middle of the plateau. Within the D expansion, results are largely 
independent of Dq. This allows us to estimate the uncertainty in v from its 
deviations from the plateau, which is the error bar quoted in the first column 
of Table |]. Expansion C also yields compatible results but, as apparent, it 
is not such a flat plateau. We find expansion B to be unreliable and the re- 
sults are not even displayed. Finally, expansion A, almost equivalent to the 
naive e expansion, shows a plateau coincident with D, but deviating slightly. 
The Hwa technique shows two extrema, one slightly above the Flory result, 
the other slightly below, so although the results are reasonable, we think we 
cannot apply it accurately as d decreases. At (2,7), the situation for the 
different extrapolations is very similar; again, expansion D gives a nice flat 
plateau consistent with the Flory estimate, and we extrapolate the best v in 
the same way as in the (2, 8) case. Expansion C gives a result completely 
consistent with D, although with not such a flat plateau, while expansion 
A starts to deviate. The cases (2, 6), (2, 5) and (2,4) follow the same trends 
as the previous ones, as apparent from Fig. [7| and the results are quoted in 
Table [I]. We conclude that expansion D is a reliable generalized e expansion 
that we can confidently apply to the physical case (2,3). Our final result is 
quoted in Table |l| Let us recall that in expansion D also gave the most 



reliable results, and it is interesting to find the same situation here. 

Concerning the ( exponent, the situation is different. We find the best re- 
sults applying Hwa's technique. For small d any of the estimates A,B,C,D ex- 
hibits a large enough plateau, and consequently, we are not confident enough 
of their robustness. The estimates quoted in Table [I] are those obtained from 
the Hwa method. As shown in Fig. [| we find two points where = 0, one 
for D < 2 the other for D > 2. Our actual estimate corresponds to the case 
Dq < 2 since, on one hand it agrees slightly better with the Flory estimate 
for large d, and on the other the curve e = seems intuitively closer to the 
actual point. 

Finally, it is reassuring that the values we obtain for v and (, although 
computed using different extrapolations, are compatible with the scaling re- 
lations Eq. |4l| 

5.2.2 Corrections to the Flory estimate 

The Flory approximation, which certainly works well for polymers, also pro- 
vides valuable insight in the case of membranes. The basic approximation 
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Figure 7: Calculation for the v exponent. The Long dashed line corresponds 
to the minimal sensitivity scheme, the dashed one is A expansion, dot-dashed 
corresponds to expansion C and the solid line is the D expansion. 



28 




29 



d 



VFlory C CFlory 



8 


0.333(5) 


0.333 


0.60 


0.600 


7 


0.374(8) 


0.375 


0.64 


0.643 


6 


0.42(1) 


0.429 


0.68 


0.692 


5 


0.47(1) 


0.500 


0.72 


0.750 


4 


0.54(2) 


0.600 


0.76 


0.818 


3 


0.62(2) 


0.750 


0.80 


0.900 



Table 1: Final results for critical exponents. 



assumes that elastic energies are comparable to self-avoiding energies. For 
self-consistency we then require that Z± = Z^. This extra condition fixes 5 
and v to be 

-4e 



4 + (D-l)(d+3) 
D + l 

VFlary ~ rf+ J • (85) 

It is also interesting to analyze the corrections to the Flory result. We can 
write at the fixed point the following equivalent definition for 5, 

6=6p - 4+(D-i)(d+3f ia log( lz ) • (86) 

Expanding to first order in e it follows that 

5 -2D 17 -AD /r ^ N 
MD) = — + jjjo^ip^. < 87 > 

where ^i(-D) is defined in Eq. ^2[ Coincidentally the extra factor appearing 
on the r.h.s of the previous equation is just 1 at D = 2, so Eq. ^ and Eq. |S7| 
give the same result at D = 2, a result already noticed in 0. 

We can apply the usual machinery to extract critical exponents. We start 
by examining the (2, 8) case. From Fig. ^| it is clear that all expansions give 
compatible results, although expansion C produces the flattest plateau. We 
use this expansion to extract our best v within this approximation. Using 
the fluctuations in the plateau, as a function of D , to estimate errors, we 
realize that results in this case are not as accurate, although still compatible 
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with previous estimates and with its Flory value. From analyzing this case 
we see also that this expansion tends to give a slight overestimate. The same 
situation holds as d decreases, as shown in Table ^. Overall, results remain 
close to the Flory estimate although with a larger uncertainty. They are still 
compatible with the values quoted in Table HJ, which we regard as our most 
accurate determinations. 



d V Vp Vy VFlory 



8 


0.333(5) 


0.34(1) 


0.34(1) 


0.333 


7 


0.374(8) 


0.39(2) 


0.39(2) 


0.375 


6 


0.42(1) 


0.44(2) 


0.44(4) 


0.429 


5 


0.47(1) 


0.51(3) 


0.51(5) 


0.500 


4 


0.54(2) 


0.60(4) 


0.60(6) 


0.600 


3 


0.62(2) 


0.71(6) 


0.70(9) 


0.750 



Table 2: Comparison of the different extrapolations for v 



5.2.3 Corrections to the Gaussian variational approximation 

Another approach that has relatively successful in dealing with problems with 
self-avoidance is the Gaussian variational approximation ||15|| . It consists in 
approximating the exact density functional by the best possible quadratic 
weight for the field h. It amounts to assuming that the self-avoiding term 
is not renormalized, that is Zj, — 1. The quantity 5 and the gyration radius 
exponent within this approximation were first computed in Q with result 

v (D-l)(d+3) 

v w = 7 j^Z§- (88) 

The value for the physical tubule is u var = |. This is clearly unphysical, being 
larger than one, but the accuracy of the Gaussian variational approximation 
should improve for large d, since it is essentially a large-c? expansion. As with 
the Flory approximation, we may determine the corrections to the Gaussian 
variational approximation within the e-expansion. From 

^- (0 -iHd + 3) 4'° e(Zt) ' (89) 
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9: Corrections to the Flory estimate for v. 
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which consistently at lowest order in e results in 



5 -2D Ut(D) 



(90) 



where z^i(-D) is defined in Eq. |82] and again the extra factor appearing on the 
r.h.s is just 1 at D = 2, so Eq. |S2| and Eq. |5U] are equal at D = 2, as reported 



approximation turns out to be the least accurate of our determinations. At 
(2, 8) all different extrapolations deliver equivalent results, but the expan- 
sion C shows the flattest plateau. The situation is the same for the cases 
(2, 7), (2, 6) and (2,5), but at (2,4) and (2,3), there are sizeable variations 
and there is no a clear plateau. For the sake of completeness, we extrapolate 
our results from this pseudo-plateau where variations are small, and quote 
an error bar from its variations, which are the results quoted in Table 0. 

6 Conclusions 

In this paper our first task was to identify the simplest free energy containing 
all the relevant operators controlling the large-distance physics of the tubular 
phase of anisotropic membranes. In this analysis essential use was made 
of rotational symmetries. Although our analysis may be modified by the 
existence of more complicated phase diagrams with non-perturbative fixed 
points in the spirit of 0, we believe that the model treated here reveals 
essential features of the physics of the anisotropic tubular phase. 

Finally we completely characterized the phase diagram and calculated 
the critical exponents by generalizing the e-expansion introduced in || . For 
the physical self-avoiding tubule we find 



in 




v 



z 



c 



0.62 



0.80 



0.75 



1.0 



-0.33 



1.65 



(91) 
(92) 
(93) 
(94) 
(95) 
(96) 
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Further improvement would necessitate a two-loop calculation for arbitrary 
D and would provide a valuable check of our extrapolation. 

These predictions may be tested via an extension of the numerical simula- 
tions described in [0] to the much more demanding model with self-avoidance. 
These simulations are currently in progress. We hope the concreteness of the 
calculations presented here will inspire further work in the rich field of the 
physics of anisotropic extended manifolds. 
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A Appendix 

In this appendix we discuss the analytical properties of several functions that 
arise in the evaluation of the quantity 1(D). We follow closely the methods 
of H , and rewrite Eq. [73] as 



with 



and 



f(u) = u 



m = t^eit /„ dzF{z) ' (97) 

F(z) = z™ r duu 2D - 3 e- zf ^ , (98) 
Jo 



4-2D 



{r+°° d , 
u J dtf-^Ki-p (t) cosf 



u 



1 r+oo D _ 3 +1/2 ~) 

+ -/ dtt—Rs-p (t) sinf )} , (99) 

2 Jo 2 u J 

where K v is a modified Bessel function. We have not been able to compute 
the integrals in Eq. |99] explicitly, except in the case D = 2. Nevertheless, 
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we know both that f(u) is a monotonically increasing function of u and its 
asymptotic behavior for small and large u. For large u we have the result, 



D 

2~ 



/(«)«-« 

while, for small u, we have 



lim/(w) 

u— »0 



- 2 rr-)rr? 



5-2D, 



l + 0(l/u)) 



(100) 



vrr 



■3-D- 



1 



lim/ (n) («) 

«— >0 



2^ r(5 - 2D) sin(f (5 - 2D)) 
, n > 



(101) 



where n stands for any derivative of u. The latter leads us to conjecture 
that the corrections to Eq. [L01| are of the type 0(e~^), as explicitly seen 
at D — 2. A plot of f(u) for different values of D is given in Fig. [TT| . 

The next step is to compute F(z) Eq. |98|. Its exact analytical form seems 
hopeless to compute, but again we can find its asymptotic limits. For small 
z we have 



F(z) 



1 V5-2D/ 



D-2 
Z5-2D 



(2 D / 2 - 2 r(l/4)r( J D/2 - 1/4))™ 5 - 2D 



[l + 0(z 



l/(5-2D)^ 



and for large z 

F(z) ~ e 
where /(0) is given by Eq. |101| . 



-f(o> 



(102) 
(103) 



D 


1(D) 


D 


1(D) 


1.6 


0.07951 


2.0 


1.263910~ 3 


1.7 


0.04643 


2.1 


1.077510" 4 


1.8 


0.02027 


2.2 


1.7155710- 6 


1.9 


0.006434 


2.3 


4.01298 KT 10 



Table 3: Sample of values for 1(D). 

The asymptotics provide valuable cross checks for the numerical integra- 
tions, and we have also used them in speeding up the numerical integration 
algorithms. A sample of values for 1(D) is provided in Table 
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Figure 11: f{u) for different values of D. The dashed line is D = 1.7, the 
dotted line D = 2.0 and the solid line D = 2.3. 
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